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NAVIGATION  WITH  HIGH- ALTITUDE  SATELLITES: 


A  STUDY  OF  RANGING  ERRORS 


I.  INTRODUCTION 

Earth  satellites  have  found  a  variety  of  useful  applications  in  a 
relatively  short  period  of  time.  One  application  for  which  satellites  seem 
uniquely  suited  is  that  of  accurate  position  determination  and  navigation. 
World-wide  navigation  using  low-altitude  satellites  has  already  been  shown 
to  be  feasible  and  an  operational  system  is  presently  being  contemplated.  1 

More  recently,  attention  has  focused  on  the  use  of  high-altitude 
satellites  for  navigation.  As  an  example,  one  possible  scheme  would  have  a 
number  of  satellites  transmitting  time  signals  which  are  generated 
independently  by  free-running  clocks  in  each  of  the  satellites.  A  user  of 
the  system  could  determine  his  position  (3  coordinates)  from  (a)  range 
difference  measurements*  for  two  pairs  of  satellites;  (b)  receiver  height 
(relative  to  the  earth’s  center);  and  (c)  the  satellite  positions.  In  this  scheme, 
the  user  is  passive,  i.  e.  ,  he  need  not  make  any  radio  transmissions  to 
navigate,  and,  furthermore,  the  user  equipment  need  not  include  a  very  stable 
clock  to  measure  the  difference  in  range  to  two  satellites.  Synchronization 
of  the  satellite  clocks  could  be  accomplished  by  having  a  control  station  on 
the  ground  which  provides  each  satellite  with  timing  corrections  relative  to 
a  master  clock  to  be  transmitted  along  with  the  actual  timing  signals  in  the 

*  Difference  in  time  of  arrival  of  signals  from  two  satellites. 
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form  of  digital  data.  Satellite  coordinates  could  also  be  included  as  additional 
digital  data  in  a  similar  manner  involving  control  stations.  A  small  number 
of  control  stations  would  be  required  to  11  see”  all  the  satellites  at  synchronous 
altitudes. 

There  are  apparent  advantages  in  using  hi gh -altitude  satellites  for 
navigation.  Position  determination  can  be  accomplished  almost  instantane¬ 
ously  and  essentially  continuously;  i.  e.  ,  the  time  required  to  obtain  a 
position  fix  or  begin  another  one  is  determined  by  the  measurement  time 
required  for  adequate  ranging  accuracy  and  the  time  to  compute  position. 

These  features  are  very  desirable  for  high-speed  users  such  as  jet  aircraft. 
Some  obvious  disadvantages  of  high-altitude  navigation  satellites  are  that 
(a)  several  satellites  are  required  for  a  position  fix  (whereas,  a  low-altitude 
system  may  require  only  a  single  satellite);  and  (b)  the  large  path  loss  from 
synchronous  altitude. 

This  report  is  concerned  only  with  the  problem  of  very  accurate 
determination  of  satellite -user  ranges  (or  range  differences)  in  a  high- 
altitude  satellite  navigation  system  such  as  that  just  described.  An  analysis 
of  navigation  errors*  using  high- altitude  satellites  has  indicated  that  10-meter 
range  accuracies  are  desirable  if  position  determinations  to  within  a  fraction 
of  a  mile  are  to  be  achieved.  We  consider  here  only  the  factors  that  would 

*  Performed  by  F.  C.  Schweppe  as  part  of  a  six-week  feasibility  study 
conducted  at  Lincoln  Laboratory  on  navigation  with  high-altitude  satellites. 

The  participants  in  the  study  were  T.  J.  Goblick,  Jr.  ,  B.  Reiffen,  F.  C. 
Schweppe,  and  R.  Teoste.  This  report  is  an  outgrowth  of  the  original  study 
of  limitations  on  range  measurements. 
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affect  ranging  systems  which  attempted  to  achieve  10-meter  accuracies, 
rather  than  the  detailed  design  of  an  optimum  ranging  system. 

There  are  two  basically  different  causes  of  errors  in  determining  the 
distance  between  a  satellite  outside  the  atmosphere  and  an  earth-based 
receiver  by  measurement  of  the  travel  time  of  radio  signals  transmitted  by 
the  satellite.  First,  errors  in  measuring  the  time  of  arrival  of  a  signal  at 
the  receiver  result  from  the  random  noise  which  masks  the  signal.  Errors 
of  this  type  depend  upon  such  parameters  as  the  radiated  signal  power, 
thermal  noise  level  due  to  the  sky  background  and  receiver  amplifier,  and 
the  shape  of  the  signaling  waveforms  used,  all  of  which  are  somewhat  under 
the  control  of  the  system  designer.  The  second  cause  of  ranging  errors  is  the 
propagation  medium  itself.  Unknown  or  randomly  varying  propagation 
anomalies  cause  changes  in  transit  time  of  radio  waves  propagating  through 
the  atmosphere.  Errors  of  this  type  may  be  minimized  by  judicious  choice 
of  operating  frequency. 

In  the  following  sections,  we  will  discuss  these  sources  of  error  in 
the  context  of  the  high -altitude  satellite  navigation  concept.  We,  therefore, 
consider  only  satellites  at  synchronous  altitude  having  earth  -  coverage 
antennas.  The  ground  receiver,  or  navigation  system  user,  must  monitor 
more  than  one  satellite  simultaneously  and  must,  therefore,  use  a  non- 
directional  antenna.  Thus  we  take  the  receiver  antenna  to  have  a  gain  of 
-3  db  (0-db  gain  with  a  3-db  polarization  loss  since  it  may  be  difficult  to 
build  a  perfectly  polarized  omnidirectional  antenna). 
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We  do  not  consider  the 
during  the  range  measurement 
of  vehicle  or  receiver  height  a 


effects  on  ranging  accuracy 
s.  We  also  do  not  consider 
s  a  ranging  problem. 


of  receiver  motion 
the  measurement 
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II.  TIME  DELAY  ESTIMATION 


Consider  the  problem  of  determining  the  time  of  arrival  or  time 

delay  at  a  receiving  station  of  a  signal  whose  shape  is  known.  If  the  signal 

transmitted  has  the  form  s(t),  the  waveform  that  the  receiver  observes 

consists  simply  of  the  transmitted  signal  delayed  by  the  time  t  that  it  takes 

for  radio  waves  to  travel  from  transmitter  to  receiver,  and  a  masking  noise, 

which  we  will  take  to  be  an  additive  white  gaussian  process  with  zero  mean 

and  ( single  -  sided)  spectral  density  N^(watts/Hz)  correlation  function 

R  (t,  u)  =  6(t-u)N  /Z~\.  The  received  waveform  from  which  t  is  to  be 
n  o 

determined  is  thus 


r  (t)  =  s  ( t  -  t  )  +  n(t )  ,  tj  <  t  Ct^,  (1) 

where  the  interval  (t^,  t^)  includes  as  much  of  the  received  data  as  is  actually 
available  for  processing.  We  have  idealized  the  problem  in  that  (a)  we  have 
not  included  any  distortion  of  the  transmitted  signal  by  the  propagating  medium; 

(b)  the  amplitude  of  the  signal  is  assumed  to  be  known  at  the  receiver;  and 

(c)  the  statistics  of  the  masking  noise  n(t)  are  completely  known.  Note  that 
r(t)  is  a  gaussian  process  with  expected  value 

E[r(t)  ]  =  s (t - t )  (2) 

and  autocovariance 

E{[r(t)  -  s(t-T)]  [r(u)  -  s(u-t)]}  =  E[n(t)n(u)]  ■  6(t-u)NQ/2  (3) 
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where  5(t-u)  is  a  unit  impulse  function  located  at  t  =  u.  Thus  the  problem  of 

determining  t  can  be  viewed  as  the  estimation  of  an  unknown  parameter  of  a 

2  3  4 

gaussian  random  process,  which  is  a  standard  statistical  problem.  1  * 

Ranging  Performance  Criterion 

Rased  on  the  observation  of  the  received  waveform  r(t),  t^  <  t  ct^* 
which  is  a  sample  function  of  a  gaussian  process,  we  must  arrive  at  an 
estimated  value  for  the  signal  delay,  denoted  t.  Thus  t  is  a  random  vari¬ 
able  dependent  on  r(t).  A  reasonable  criterion  by  which  to  judge  various 
estimation  schemes  is  that  of  the  mean  squared  error  between  the  estimate 
t  and  the  actual  value  of  t;  i.  e.  , 

eT2  =  E(t  -  -)2  .  (4) 

If  we  have 

E(t)  =  t 

*  2 
then  t  is  called  an  unbiased  estimator.  In  this  case,  is  just  the 

variance  of  the  estimate  t. 

Bound  on  Estimation  Accuracy 

Next  we  will  find  the  minimum  variance  of  any  estimate  for  t  (the 

4 

Cramer-Rao  lower  bound  on  the  variance  of  an  estimate  ).  In  order  to 
write  a  probability  density  for  the  received  waveform  r(t),  we  first  expand 
r(t)  in  the  Fourier  series 
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(5) 


where 


00 

r(t)  =  7  ri  cp.  (t)  ,  tj  <  t  <t2  , 

i  =  1 


r . 
i 


r (t)  OX  (t)  dt 


a  nd 


Using  Eq.  1,  we  can  write 


where 


a  nd 


r.  =  s .  (  t  )  +  n. 
ii  i 


s .  (  t  ) 


s(t- t)  cry  (t)  dt 


n.  = 

l 


n(t)  cp.  (t)  dt 
i 


(6) 


(7) 


(8) 


(9) 


(10) 


If  s(t)  is  restricted  in  bandwidth  by  requiring  that  only  a  small  fraction  of  its 
energy  lies  outside  a  specified  range  of  frequencies,  many  of  the  coefficients 
s^(t)  will  be  restricted  to  be  very  small.  This  implies  that  s(t-^)  can  be 
represented  by  a  truncated  Fourier  series 

M 

s(t--r)  =  Yj  S.(t)  cp.  (t)  ,  t  J  C t;  C 1 2 .  (11) 

i  =  l 
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without  significant  error.  Since  n(t)  was  taken  to  be  white  gaussian  noise, 
the  coefficients  m  are  statistically  independent  regardless  of  the  set  of  ortho¬ 
normal  functions  fcp^  (t)]  that  are  used.  Since  s^(t)  ~  0  for  i  >  M,  the 
coefficients  m  for  i  >  M  are  irrelevant  to  the  problem  of  estimating  t  since 
they  do  not  yield  any  information  about  any  coefficients  for  i  <  M.  Thus  we 
may  safely  truncate  the  expansion  (5)  as 


M 

r(t)  =  V  r  CD.  (t)  ,  t  <t  <t.  .  (12) 

.* — 11  i  4 

1=1 


Thus  r(t)  is  specified  by  the  set  of  Fourier  coefficients 


r  =  (r  ,  r ,  i*M)  (  1  3) 


for  which  we  can  write  a  probability  density  conditional  on  t  as 


p(r  |  t) 


n 


i*l 


expf  -  ^r .  -  s. 

li 


(hNq) 


(t)1z/No} 

T/2— 


(14) 


Denoting  our  estimate  as  t  (r)  since  it  is  a  function  of  the  observables 

r.,  the  bias  of  the  estimate  is  defined  as 
i  - 


'Ht)  =  J  [t  (r)  —  T  ]  p(r  |  t)  dr 


(15) 


Differentiating  with  respect  to  t  ,  we  get 


djrM 

d  t 


=  i'  (r) 


Sp(£  |  t) 

J  r t  <£.)  *  T  3  — ^ -  dr  -  ] . 


(16) 
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U  sing 


d  lnp(r|^) 

we  get  from  the  Schwarz  inequality 


j  Op(r|T) 
p(r  |t)  d  t 


d  In  p(r  |  t) 

1  +  V  (t)  =  J  [T(r)-r]  - -  P (£  |  t )  dr 

_  d  In  p(r  |t) 

=  J  {rT(r)-T]  /p(r  |  t)}  f g— -  / p(r  |  t)}  dr 


f  j  ['(r)-T]2  p(r  |t)  dr }"{  f[ 


Slnp(r 


|t) 


p  ( r  |  t  )  d  r  }  ' 


and  hence 


2 


e 

T 


[i  +V  (th 

5  In  p(r  I 


Ef 


5  T 


2 

tT 


2 

] 


(17) 


This  is  a  lower  bound  to  the  mean  squared  error  for  any  kind  of  estimate  of 
t  and  is  satisfied  with  equality  if 


5  In  p(r  |  r) 

- 5— -  =  M^)  [r(r)-T] 


where  1c(t)  is  independent  of  t. 


(18) 
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Bound  on  Range  Accuracy 


A  lower  bound  on  the  mean  squared  error  with  which  the  time  of 
arrival  of  a  signal  can  be  estimated  by  any  scheme  can  be  found.  First  note 
that  from  Eqs.  7,  11,  and  12, 


t 


j 


z 

1 


r  (t)  —  s(t-  t)  dt 


(19) 


Using  this  relation  in  Eq.  14  results  in 


^  In  n(r 


[  r(t)  -  s(t-T)l 


ds(t-T) 

3“t 


dt . 


(20) 


Sine  e 


8s(t-  t)  _  ds(t-  t) 

St  St 


(21) 


we  use  this  together  with  Eqs.  1  and  3  to  show  that 


3  lnp(r  |t)  ,2 

Br  J 


N 


dt 


B  s(t- t) 

“St 


J  du  5s5(^~t)  E ( n(t)  n(u)] 


2 

"  TT- 

O 


r» 


I  Bs(t-T) 

I  Tt 


(22) 
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The  Fourier  transform  of  the  signal  s(t-T)  is 


S(f) 


s ( t - t )  e~^  ^  dt  , 


(f  in  Hz) 


(23) 


so  that  the  transform  of  the  signal 


Ss(t-T) 

~5t - 


is  obtained  by  integrating  by  parts, 


and  is 


t2 

J  d-S^~'r-  e'}2rTit  dt  =  j 2tt f  S(f).  ( 24) 

*1 


From  Parseval’s  theorem  the  received  signal  energy  can  be  written  as 


E 

s 


■  J  |S(f)|2df 


(25) 


and  similarly, 


j 

t 

i 


(  2tt)  2  J  f2  |S(f)  |  2  dt  =  (ZttB)ZE  (26) 
-00 


where 


J  f2  |  S(f)  |  2  df 
-  00 

J  |S(f)|2df 


(27) 


can  be  interpreted  as  the  mean  squared  bandwidth  (in  Hz  )  of  the  (normalized) 
spectrum  of  the  ranging  signal  s(t).  Thus  the  mean  squared  ranging  error 
can  be  lower  bounded  as 


2 


[1  +  *  (t)1 

IE 


N 


-  (ZttB)2 


(28) 


The  parameters  of  interest  regarding  ranging  accuracy  appear  to  be 
the  receiver  signal -to -noi  se  ratio  E  /N ^  as  well  as  the  (mean  squared)  band¬ 
width  of  the  ranging  signal.  However,  it  is  not  yet  apparent  that  Eq.  Z8 
represents  a  useful  bound  on  ranging  accuracy. 

Maximum  Likelihood  Estimators 

It  is  impossible  to  derive  the  best  estimator  for  t  in  the  mean  squared 
error  sense  in  the  general  case.  Instead,  the  maximum  likelihood  estimate 
(MLE)  for  t  (in  the  absence  of  an  a  priori  probability  density  for  t)  is 
commonly  substituted  with  good  reason  (as  will  presently  be  shown);  i.  e.  , 
that  value  of  t  that  maximizes  p(t|t)  is  chosen  as  the  estimate  of  t.  Thus 
the  MLE  satisfies 


S  In  p( r  |t) 


=  0 


a  t 


and  using  Eq.  19,  we  obtain  an  integral  equation  for  t  as 


j  ( r(t)  -  s(t-T)]  dS-(*-T)-  dt  =  0 


a  t 


(29) 


(30) 
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Note  that  in  the  case  of  no  noise  mashing  r(t),  r(t)  =  s(t-r)  and  t=  t 
satisfies  Eq.  30.  In  the  case  of  very  low  noise  power  in  the  signal  bandwidth, 
p(t|t)  will  have  a  narrow  peak  in  the  vicinity  of  T.  If  the  width  of  this  peak 
is  small  enough,  we  may  approximate  s(t-T)  by  the  expansion 


S(t-T)  «  S(t-T)  +  (T-T) 


as(t-T) 


(31) 


allowing  Eq.  1  to  be  rewritten 


r (t)  —  s (t  —  t )  «  (t-t)  — —  +  n(t ) . 


Using  Eqs.  26  and  32  in  30  yields 


(32) 


T  + 


( 2nB)u  E  t, 

s  1 


S  nw 


ds(t-  t) 

3  t 


dt. 


(33) 


Since  n(t)  was  assumed  to  have  zero  mean, 


E(t) 


(34) 


and  the  MLE  is  unbiased  for  all  t  in  this  so-called  weak  noise  case.  There¬ 


fore,  i!i  (t)  <3-  0.  The  variance  of  the  estimate  is  thus 
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2  r/*  x2 

e  =  E(t-t) 

T 


- '  /2dtMid)  JZ  du55(Hll)E[n(t)n(u)l 

U^BpE^  tj  Sr  tj  St 


TE 


(35) 


N 


-  (2tB)2 


which  shows  that  in  this  weak  noise  case  the  MLE  of  t  is  efficient,  in  the 


sense  that  it  achieves  as  small  a  variance  as  any  estimate  for  t. 
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III.  IDEAL  RANGING  PERFORMANCE 


Ideal  Ranging  Signal 

If  the  ranging  signal  is  restricted  to  the  band  of  frequencies 
f  <  W(in  Hz)* 

then  the  problem  of  finding  the  signal  with  maximum  mean  squared  bandwidth 
(with  normalized  spectrum)  is  analogous  to  that  of  finding  the  distribution  of 
a  given  amount  of  mass  along  a  finite  length  rod  to  achieve  the  maximum 
moment  of  inertia.  The  optimum  signal  in  this  sense  would  have  a  spectrum 
S(f)  whose  area  was  concentrated  at  the  edges  of  the  allowable  band  and  w'as 
zero  elsewhere  within  the  band,  the  resulting  mean  squared  banduddth  being 
of  course.  This  is  approximately  achieved  by  a  sinusoidal  signal 
of  frequency  W  and  duration  T  seconds  in  which 

w  »  1  . 

The  sinusoid  can  thus  be  considered  an  ideal  ranging  signal  in  the  sense  that 
it  has  maximum  mean  squared  bandwidth  for  a  given  spectral  occupancy. 
According  to  Eq.  35,  a  simple  radio  frequency  (RF)  carrier 

s  (t)  =  A  co  s  (2tt  f  t  +  f6  ) 

o  c 

can  be  located  in  time  with  great  accuracy  (if  the  parameters  f  and  $  are 

o  o 

accurately  known  at  the  receiver).  It  is  not  even  necessary  to  do  the  range 

*  All  frequencies  and  bandwidths  are  expressed  in  Hz  (cycles  per  second). 
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measurement  at  the  RF  frequency  f  since  multiplication  of  the  received 
signal  by  Z  cos  Zrr(f^  —  f^)t  gives 

Zs(t-T)cos  Ztt  ( f  —  f  )t  =  Af  cos(Znf  t  —  Ztt  f  t  +  0  ) 
o  1  ]  o  o 

-I-  COs[Zrr(Zf  — f,)t  — Zrrf  t  +  0  )] 
O  J  o  o 


and  bandpass  filtering  can  then  eliminate  the  term  at  frequency  Z f  f  ^ .  The 
result  of  this  transformation  of  the  signal  is 

f 

A  cos  [2nf  (t  —  •£—  t)  +  . 


Calling  f  r/f  =  ti,  the  accuracy  of  the  estimate  of  T|  is  inversely  proportional 


to  f^  ,  but 


,  f  Z  7 
Z  /  o  »  Z 

en  *  V  ^ 


and  so  is  really  independent  of  f ^ *  This  is  a  satisfying  situation  since 

it  indicates  that  the  range  accuracy  is  inherently  related  only  to  the  carrier 
f r  equency . 

Resolving  Ambiguities 

There  is  another  consideration  with  regard  to  ranging  signals  that 
we  have  not  yet  touched  upon.  That  is  the  question  of  range  ambiguities. 
Observation  of  a  CW  sine  wave 


sin  Ztt  W  ( t  -  t  )  =  sin(2TTWt-^) 


*  This  was  pointed  out  by  B.  Reiffen. 
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cannot  be  used  to  determine  the  phase  angle  0  uniquely.  There  will  always 
be  an  ambiguity  of  some  multiple  of  Zn  radians  in  the  determination  of 
implying  that  the  time  delay  r  will  always  be  ambiguous  by  multiples  of  the 
period  of  the  sine  wave.  This  effect  implies  that  there  are  multiple  solutions 
to  Eq.  30,  a  question  that  we  did  not  consider  before. 

A  sine  wave  is  still  an  efficient  signal  with  which  to  achieve  high 
range  accuracy  as  long  as  a  separate  scheme  is  first  used  to  reduce  range 
uncertainty  to  less  than  one  period  of  the  sine  wave,  thus  resolving  the 
ambiguity  problem  by  ruling  out  all  but  one  of  the  many  solutions  of  Eq.  30. 
Many  practical  ranging  systems  use  this  general  scheme  of  a  separate  signal 
for  coarse  range  determination  to  resolve  ambiguities  and  a  sine  wave  for  the 
high  resolution  or  fine  range  measurement.  The  range  accuracy  of  such 
hybrid  systems  can  still  be  computed  from  Eq.  35  by  taking  B  and  Eg  to  be 
the  frequency  and  the  energy  of  the  fine-range  sine  wave  respectively. 

RF  Bandwidth  Occupancy  (Coarse-Fine  Measurement  Systems) 

We  face  a  dilemma  in  trying  to  relate  range  accuracy  to  RF  bandwidth 
occupancy  of  the  ranging  signal  since  we  have  already  noted  that  an  unmodu¬ 
lated  RF  carrier  can  be  used  for  the  high  resolution  range  measurement, 
which  occupies  essentially  no  bandwidth.  In  order  to  remove  range  ambigu¬ 
ities,  a  coarse  range  signal  will  have  to  be  sent  (at  RF)  and  so  there  will  be 
some  bandwidth  occupancy.  Suppose  the  fine  range  measurement  is  made  using 
the  carrier  of  frequency  fQ.  The  coarse  range  signal  need  only  provide  an 
estimate  of  time  of  arrival  with  enough  accuracy  to  resolve  the  fine  range 
ambiguities.  The  root  mean  squared  (RMS)  error  of  the  coarse  measurement, 
denoted  a  ,  can  thus  be  taken  as 
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(36) 


o  -  j-  ,  0  <  a  <  1  • 

c  r  1 

o 

Therefore,  for  a  given  signal  energy,  the  mean  squared  bandwidth  of  the  coarse 
range  signal  should  be  (using  Eq.  36  in  Eq.  35), 

,  f  2 

B  =  -  (37) 


cr 


ZE 

T9 


cr  (Zn  a)2 


while  the  system  accuracy  is  given  by 


ZE 

N  o 


(38) 


where  E^  and  are  the  carrier  and  coarse  range  signal  energies  respec¬ 


tively.  Combining  these  equations  gives 


B 


1 


c r  7  ZE ZE 

o  cr  Z  Z 

(2n1  "N - W~  a  eT 

o  o 


(39) 


In  principle,  we  can  find  no  objection  to  having  the  same  signal  serve  as 
both  the  coarse  and  fine  ranging  signal.  The  fine  range  measurement  could  be 
made  using  the  fine  structure  or  carrier  of  the  signal  at  RE  or  after  transla¬ 
tion  to  an  intermediate  frequency  (IE).  If  the  translated  signal  is  narrowband, 
we  can  take 

B2  Z 

o 

for  the  fine  measurement.  The  coarse  measurement  would  effectively  use  the 
envelope  or  carrier  modulation.  If  this  modulation  has  mean  squared  bandwidth 
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B  ,  then  the  RF  signal  can  occupy  an  RF  bandwidth  of  no  less  than  B 
cr  °  7  cr 

[and  this  bandwidth  can  be  realized  only  if  (a)  the  modulation  is  transmitted 
single- sideband,  and  (b)  the  actual  occupied  bandwidth  of  the  modulating  wave¬ 
form  is  essentially  the  minimum  possible  value]  .  In  this  case,  we  would 

have 

E  =  E  =  E  .  (40) 

o  c  r  s 


Since  we  must  have 

E  <  E 
o  s 

E  <  E 
c  r  s 


(40) 


use  of  Eq.  40  in  Eq.  39  results  in  a  formal  lower  bound  on  the  required  RF 
bandwidth  for  any  ranging  system  as 


W  £ 
RF 


1 

,  >2  s 

(2n)  -jj-  acl 


(41) 


The  accuracy  of  the  fine  measurement  would  still  be  given  by  Eq.  38  (with 

E  =  E  ). 
o  s 

The  above  discussion  is  intended  merely  to  point  out  that  (a)  the  desired 
range  accuracy  alone  does  not  imply  anything  about  the  occupied  bandwidth 
W  of  the  ranging  signal;  (b)  the  a  priori  range  uncertainty  is  related  to  the 
required  RF  bandwidth  of  the  ranging  signal;  (c)  a  lower  bound  to  was 

given  for  any  ranging  scheme  with  large  enough  a  priori  range  uncertainty  to 
necessitate  a  coarse  range  measurement.  Study  of  some  examples  shows  that 
for  a  given  si gnal-to-noi se  ratio,  the  lower  bound  on  of  Eq.  41  can  be 

approached  well  within  a  factor  of  10  in  certain  situations. 
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that  satisfy  Eq.  39 


In  Table  I  some  combinations  of  E  /N  and  W 

s  o 

-  8 

with  equality  have  been  computed  for  =  10  sec  (3-meter  RMS  range 
error).  The  parameter  a  essentially  controls  the  probability  that  the  coarse 

measurement  resolves  the  fine  range  ambiguities  incorrectly,  i.  e.  ,  the 
probability  of  a  grossly  inaccurate  or  anomalous  range  determination.  Reason¬ 
able  values  for  a  would  be  between  0.  5  and  0.  1,  and  for  purposes  of  illustration 
we  will  choose  a  =  0.  3  (-5  db). 


Table  I 


E 

<db> 

W~  in  KHz 

K  x 

o 

6 

100 

16 

10 

26 

1 

36 

0.  1 

Double-Sideband  (DSB)  Amplitude  Modulation 

By  way  of  illustrating  a  more  conservative  approach  to  obtaining  the 
power-bandwidth  requirements  for  ranging,  we  next  discuss  systems  in  which 
the  ranging  signal  sQ(t)  is  transmitted  in  the  form 

s  (t)  =  s  (t)  cos  2tt  f  t 

If  the  carrier  parameters  are  all  known  and  removal  of  the  carrier  (demodula 
tion)  can  be  accomplished  without  degrading  the  modulation  (for  ranging  pur¬ 
poses),  the  range  estimation  can  be  done  using  the  waveform  SQ(t).  This  can 
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be  accomplished  in  practice  for  a  variety  of  modulations  (for  instance,  those 
which  allow  carrier  removal  with  a  phase-locked  loop).  Thus  if  sQ(t)  has 
mean  squared  bandwidth  B^,  then  it  must  occupy  a  low-pass  bandwidth  of  at 
least  BHz,  which  implies  that  the  bandwidth  occupied  by  s(t)  satisfies 


W 


DSB 


:>  ZB. 


Using  this  in  Eq.  35  results  in  a  bound  on  performance  of  systems  which  use 
only  the  double  -  sideband  modulation  for  ranging.  An  example  of  such  systems 
is  a  short-pulse  radar  system  in  which  the  pulse  width  governs  the  ranging 
accuracy.  The  bound  on  performance  is 


W 


z 

DSB 


o 


(4Z) 


and  this  bound  is  plotted  in  Fig.  1  (again  for  =  1 0  sec).  Denoting  the 
average  power  of  the  received  signal  as 


P 

r 


(43) 


we  can  write 


E  P 

s  _  r 

N - 

o  o 


which  brings  out  explicitly  the  measurement  time  T  as  a  parameter.  In 
Fig.  1  the  Eg/N  axis  is  thus  also  labeled  as  P^/N^  corresponding  to  a 
measurement  time  of  1  sec. 


Z1 


Es/N0(dB)  AND  Pr/N0(dB)  FOR  T  =  \  sec 


Fig.  1.  Ideal  DSB-AM  ranging  system  trading  curve  for  bandwidth  vs 
signal-to-noise  ratio. 


22 


Note  the  radical  difference  in  the  required  si gnal - to - noi se  ratio  -  band  width 


combinations  between  the  ideal  DSB  system  and  those  given  in  Table  I  for  the 

same  range  accuracy.  The  DSB  requirements  are  conservative  in  technology, 

although  they  still  correspond  to  ideal  performance  in  the  sense  that  we  took 

W  =2B  The  bound  on  occupied  bandwidth  in  a  DSB- AM  ranging  system  can 
DSB  * 

be  approached  well  within  a  factor  of  10,  for  a  given  signal -to  -  noi  se  ratio  and 
a  large  a  priori  range  uncertainty.  Even  this  margin  does  not  lead  to  prohibi¬ 
tive  satellite  power  levels  and  bandwidths.  The  performance  indicated  in 
Table  1  may  perhaps  be  approached  as  closely  with  further  work  on  the  type 
of  signal  processing  required. 


Satellite  Power  Requirements 

E  P  T 
s  r 

Having  calculated  the  minimum  required  si  gnal  -  to  -  noi  se  ratio  =  N - 

o  o 

at  the  receiver  corresponding  to  a  given  signal  bandwidth  for  a  DSB-AM 
ranging  system,  we  now  translate  this  into  satellite  transmitter  power  require¬ 
ments.  As  an  example  calculation,  the  signal  power  P^  received  by  an 
antenna  with  gain  =  -3dB  from  a  satellite  at  synchronous  altitude  with 

an  earth  coverage  antenna  (gain  G  =  18  dB)  at  an  operating  frequency  of 

sat 

300  MHz  is 


P  (dBw) 
r 


Psat  Gsat 


path 


G 

rec 


=  P  (dBw)  +  18  dB  -  174  dB  -  3  dB 

S  St  L 


=  P  „  (dBw)  -  159  dB 
sat 


(44) 
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where  P  (dBw)  is  the  satellite  radiated  power  expressed  in  decibels  re 
sat 

1  watt.  The  path  gain  varies  inversely  as  the  square  of  the  operating  frequency* 
so  can  be  regarded  as  a  function  of  frequency  for  a  fixed  Pg  . 

The  noise  power  density  is  due  to  sky  background  noise  and  receiver 

noise,  both  of  which  are  dependent  on  operating  frequency.  The  noise  tempera¬ 
ture  of  cosmic  noise  in  the  direction  of  the  galactic  plane  is  plotted  as  a  func¬ 
tion  of  frequency  in  Fig.  2.  Also  shown  in  this  figure  are  the  noise  temperature 

5 

of  an  uncooled  parametric  front  end  and  the  sum  of  the  sky  and  parametric 

receiver  noise  temperatures.  The  right-hand  scale  of  Fig.  2  translates  these 

noise  temperatures  to  noise  power  densities. 

Returning  to  our  example  calculation  of  P  ,  the  receiver  noise 

r  sat 

power  density  is  seen  from  Fig.  2  to  be 


N  =  -200  dBw  (re  1  Hz) 


at  the  assumed  operating  frequency  of  300  MHz.  If  an  RF  bandwidth  occupancy 
5 

of  10  Hz  is  chosen,  Fig.  1  indicates  that  we  must  have 

P 

r  i<7dB 

o 

corresponding  to  a  measurement  time  T  =  ]  second.  Therefore,  using  Eq.  44, 

P 

^  (dB)  =  P  (dBw)  -  N  (dBw)  —  1 59  dB 
N  sat  o 

o 

and  so 

P  ^  47  dB  —  200  dBw  +  159  dB  =  6  dBw  (4  watts), 
sat 
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NOISE  TEMPERATURE  (°K) 


-190 


-195 


-200 


-  205 


_ L _ _ I _ I _ I  1  t  1  J _ I _ I _ I _ I _ 1111 

100  200  400  1000  2000  4000 

OPERATING  FREOUENCY  (MHz) 

Fig.  2.  Thermal  noise  environment  of  ranging  system. 
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NOISE  POWER  DENSITY  N  d B w  (re  1  Hz) 


Four  watts  of  radiated  satellite  power  are  required  to  achieve  an  RMS  ranging 

error  of  3  meters  at  a n  operating  frequency  of  300  MHz.  At  a  frequency  of 

1Z00  MHz,  the  path  gain  decreases  by  1Z  dB  and  is  about  -Z08  dBw,  a 

decrease  of  8  dB  (see  Fig.  Z),  implying  that  P  >  6  dB  +  1Z  dB  —  8  dB  =  10  dBw 

sat 

(10  watts).  By  scaling  path  gain  in  this  manner  and  using  Fig.  Z,  the  required 

P  can  be  obtained  for  all  operating  frequencies  of  interest  corresponding  to 
sat 

=  47 dB 

o 

This  curve  has  been  plotted  in  Fig.  3.  The  fact  that  this  value  of  P^/N^ 

5 

corresponds  to  =  10  Hz  in  an  ideal  DSB  ranging  system  is  also  indi¬ 

cated  on  the  curve.  By  using  Fig.  1  to  obtain  the  required  P^/N^  for  other 

values  of  the  ideal  system  bandwidth,  curves  of  P  .  can  be  calculated  for 

7  sat 

various  ideal  DSB  ranging  system  combinations  of  E  /N  and  Such 

”  &  7  s  o  DSB 

curves  also  appear  in  Fig.  3  and  it  should  be  noted  that  they  apply  to  any 
(non-ideal)  DSB  ranging  system  that  requires  the  indicated  P  /N^  (but  perhaps 
a  larger,  non-ideal  bandwidth). 

For  a  given  P^/N  ,  a  minimum  of  satellite  power  is  required  for 
operating  frequencies  in  the  vicinity  of  300  MHz,  but  this  is  a  very  broad 
minimum.  The  curves  rise  for  low  frequencies  because  the  cosmic  noise 
increase  requires  greater  satellite  power.  The  curves  rise  again  for  high 
frequencies  because  of  the  increasing  and  the  decreasing  aperture  of  the 

receiving  antenna  with  increasing  frequency.  However,  operating  frequencies 
well  above  300  MHz  are  actually  desirable  in  order  to  avoid  frequency- 
dependent  ionospheric  propagation  anomalies. 
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SATELLITE  TRANSMITTER  POWER  (dBw) 
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Fig.  3.  Satellite  transmitter  power  required  for  "Ideal"  DSB 
ranging  system  vs  operating  frequency. 

Assumptions: 

1  .  Synchronous  altitude  satellite 
with  earth  coverage  antenna 

2.  Parametric  receiver  front  end 

3  .  Receiver  antenna  gain  —  3  dB 

4.  Measurement  time  1  second 

5.  RMS  range  error  3  meters 
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IV.  ATMOSPHERIC  PROPAGATION  ANOMALIES 


Propagation  anomalies  introduced  by  the  atmosphere  can  be  corrected 
out  to  the  extent  that  they  are  predictable.  Of  course,  random  anomalies 
cannot  be  so  removed  and  thus  represent  fundamental  limitations  imposed  on 
ground-based  receivers  that  attempt  to  make  very  precise  measurements  of 
signal  parameters. 

The  Ionosphere 

The  inhomogeneities  of  the  ionosphere  result  in  a  refractive 
bending  of  radio  waves  which  is  frequency  dependent.  This  bending  of  the 
propagation  path  results  in  a  difference  between  the  direction  of  arrival  of 
radio  waves  at  the  receiver  and  the  actual  direction  of  the  transmitter.  This 
effect  does  not  concern  us.  Ionospheric  refraction  also  results  in  a  propaga¬ 
tion  path  length  which  is  different  from  the  actual  distance  between  trans¬ 
mitter  and  receiver.  This  effect  is  the  one  with  which  we  are  concerned  in 
measuring  time  of  arrival  of  signals  from  exo-atmospheric  satellites. 

Millman^  has  computed  the  range  error  due  to  ionospheric  refraction 
for  a  model  ionosphere.  His  results  are  shown  in  Fig.  4  for  an  exo-atmos¬ 
pheric  transmitter  and  an  earth-based  receiver.  Note  the  difference  in  range 
errors  due  to  the  different  elevation  angles  (angle  above  the  horizon)  of  the 
transmitter.  This  variation  with  elevation  angle  is  swamped  by  the  day-night 
variation  of  the  electron  density  in  the  ionosphere,  which  would  be  difficult 
to  predict.  This  data  indicates  that  frequencies  above  1000  MHz  would  allow 
overall  range  accuracies  of  the  order  of  10  meters  RMS  even  if  ionospheric 
refraction  is  ignored. 
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IONOSPHERIC  RANGE  ERROR  (ft) 


Fig.  4.  Range  errors  due  to  ionospheric  refraction  for  a  one-way 
path  from  a  high-altitude  satellite  to  the  earth's  surface. 
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The  Troposphere 


Millman  has  also  calculated  the  range  errors  due  to  tropospheric 
refraction  for  a  standard  atmosphere  for  conditions  of  0#  and  100#  humidity 
(see  Fig.  5).  These  errors  depend  on  air  temperature  and  humidity  but  are 
independent  of  operating  frequency.  If  corrections  for  humidity  and  tempera¬ 
ture  are  not  made,  the  errors  due  to  these  causes  could  be  considered  as 
random  ones.  For  elevation  angles  below  4°,  the  errors  due  to  neglecting 
humidity  are  greater  than  10  feet.  Additional  errors  due  to  neglecting  tem¬ 
perature  are  of  the  order  of  15  to  20#  of  the  humidity  errors.  Thus,  if  simple 
range  corrections  are  made  for  tropospheric  refraction  (ignoring  temperature 
and  humidity  corrections),  the  elevation  angles  at  which  satellites  are  usable 
for  navigation  (i.  e.  ,  resulting  in  range  errors  of  the  order  of  10  meters  RMS) 
is  restricted  to  be  above  8°. 

Millman  also  points  out  that  refraction  is  a  function  of  air  density  and, 
hence,  most  of  the  range  error  is  due  to  radio  path  refraction  at  low  altitude. 
High-flying  aircraft,  therefore,  would  have  to  modify  their  t ropo  -  r ef raction 
correction  or,  perhaps,  might  even  be  able  to  ignore  it. 
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TROPOSPHERIC  RANGE  ERRORS  (ft) 
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3-66-5485 


Fig.  5.  Range  errors  due  to  tropospheric  refraction  for  a  one-way 
patn  from  a  high-altitude  satellite  to  the  earth’s  surface. 
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V.  MULTIPATH 


A  receiving  antenna  with  a  wide-angle  pattern,  perhaps  a  hemispherical 
pattern,  is  needed  to  receive  the  timing  signals  from  three  satellites  simul¬ 
taneously,  especially  since  the  satellites  cannot  be  clumped  in  one  part  of  the 
sky  for  accurate  position  determination.  Such  an  antenna  installed  on  an  air¬ 
craft  may  also  receive  the  satellite  signals  reflected  from  the  earth  (multi- 
path  propagation)  with  deleterious  effects  on  ranging  accuracy.  Severe  multi- 

path  propagation  of  satellite  signals  has  been  observed  at  an  operating  fre- 

9 

quency  of  229  MHz  in  aircraft  flying  over  water  and  smooth  terrain.  Studies 
are  presently  being  conducted  to  determine  its  effect  on  digital  communication 
links  to  aircraft  via  satellites.  The  much  more  precise  measurements  required 
of  a  satellite  navigation  receiver  will  be  more  seriously  affected  than  a  digital 
communication  link.  Hence,  more  refined  experiments  will  be  needed  to 
determine  the  effects  of  multipath  on  ranging  accuracy. 

The  angle  between  the  waves  arriving  at  the  receiver  via  the  direct 
path  and  the  ea rth  -  reflected  paths  decreases  as  the  satellite  elevation  angle 
decreases.  Also  the  reflections  are  expected  to  be  strongest  for  low- 
incidence  angles.  An  ideal  up-looking  antenna  with  a  hemispherical  pattern 
would  not  be  sensitive  to  signals  other  than  those  arriving  directly  from  a 
satellite.  A  practical  antenna  would  have  to  have  a  sharp-edged  pattern  to 
attenuate  the  multipath  signals  well  for  satellites  at  low  elevation  angles. 

This  may  be  difficult  to  achieve  in  a  practical  antenna,  especially  in  the  UHF 
band.  Multipath  considerations  may,  therefore,  also  restrict  the  use  of 
satellites  at  low  elevation  angles. 
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At  higher  frequencies  (having  wavelengths  of  the  order  of  a  foot  or  less), 
even  an  ocean  surface  would  appear  to  be  a  "rough  scatterer"  of  radio  waves. 

One  would  expect  that  reflected  signals  will  be  incoherent,  thus  simply  contrib¬ 
uting  to  the  receiver  noise  level  masking  the  direct  path  signals.  The  effect  of 
incoherent  multipath  would  be  much  less  serious  than  specular  or  coherent 
reflections.  The  antenna  design  problems  may  also  be  more  readily  solved  at 
the  higher  frequencies.  Multipath  considerations,  therefore,  seem  to  favor  the 
use  of  operating  frequencies  above  the  UHF  band  for  several  reasons. 

The  effects  of  multipath  on  ranging  signals  can  be  experimentally  studied 
rather  conveniently  if  a  directional  antenna  is  available  to  provide  a  multipath- 
free  reference  signal  from  which  the  effects  of  satellite  motion  and  satellite 
clock  instability  can  be  determined.  Such  a  reference  signal  can  then  be  used 
with  a  down-looking  antenna  to  study  characteristics  of  the  multipath,  such  as 
strength,  coherence,  and  bandwidth.  The  effects  of  multipath  on  ranging 
accuracy  could  also  be  directly  determined  for  an  experimental  receiving 
antenna  with  such  a  reference  signal  by  simply  comparing  the  ranges  measured 
at  the  two  antennas. 

Since  many  high -accuracy  ranging  systems  do  the  fine  range  measure¬ 
ment  on  a  CW  sine  wave,  any  reasonably  stable  satellite  oscillator  could  pro¬ 
vide  some  useful  data  with  the  combination  of  receiving  antennas  described 
above.  Navy  navigation  satellites  could  provide  some  data  at  150  and  400  MHz. 
An  experimental  satellite  transmitting  timing  signals  from  a  very  stable  clock  is 
currently  being  contemplated.  Such  a  satellite  should  be  able  to  provide  (with 
suitable  ground  equipment)  almost  all  the  necessary  information  at  its  operating 


33 


frequency  (probably  in  the  UHF  band).  There  is  a  need  for  measurements  in 
the  800-  to  2000-MHz  range,  especially  since  this  appears  to  be  the  most 
promising  one  for  navigation.  If  a  satellite-borne  clock  is  not  available  in 
this  frequency  range,  signals  transponded  by  a  satellite  may  suffice.  Such 
features  as  switched  antennas  in  the  satellite  must  be  avoided  however. 
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VI.  SUMMARY 


1.  It  appears  that  the  main  uncertainty  concerning  the  accurate  measure¬ 
ment  of  range  to  high-altitude  satellites  is  the  effect  of  multipath  reception. 
Experimental  investigations  are  necessary  to  resolve  the  issue  before  a  firm 
statement  can  be  made  with  regard  to  the  limitations  on  range  accuracy  for  a 
high  -  altitude  satellite  navigation  system. 

Z.  Ionospheric  and  tropospheric  refraction  do  not  seriously  affect 
range  accuracy  for  operating  frequencies  above  1000  MHz  and  satellites  at 
least  8°  above  the  horizon  (Figs.  4  and  5). 

3.  Satellite  transmitter  power  levels  and  bandwidths  consistent  with 
a  range  accuracy  of  3  meters  RMS  are  not  prohibitive  (Fig.  3). 

4.  Further  work  on  ranging  signal  design  is  needed,  taking  into 
account  the  complexity  of  satellite  and  user  equipment  and  sensitivity  to 
multipath,  as  well  as  satellite  power  and  bandwidth  requirements. 
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